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For a family of logistic- like maps, we investigate the rate of convergence to 
the critical attractor when an ensemble of initial conditions is uniformly spread 
over the entire phase space. We found that the phase space volume occupied 
by the ensemble W{t) depicts a power-law decay with log-periodic oscillations 
reflecting the multifractal character of the critical attractor. We explore the 
parametric dependence of the power-law exponent and the amplitude of the 
log-periodic oscillations with the attractor's fractal dimension governed by 
the inflexion of the map near its extremal point. Further, we investigate the 
temporal evolution of W{t) for the circle map whose critical attractor is dense. 
In this case, we found W{t) to exhibit a rich pattern with a slow logarithmic 
decay of the lower bounds. These results are discussed in the context of 
non-extensive Tsallis entropies. 
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I. INTRODUCTION 



Nonlinear low-dimensional dissipative maps can describe a great variety of systems with 
few degrees of freedom. The underlying nonlinearity can induce the system to exhibit a 
complex behavior with quite structured paths in the phase space. The sensitivity to initial 
conditions is a relevant aspect associated to the structure of the dynamical attractor. In 
general, the sensitivity is measured as the effect of any uncertainty on the system's variables. 
For systems exhibiting periodic or chaotic orbits, the effect of any uncertainty on initial 
conditions depicts an exponential temporal evolution with ^(t) = limAx(o)^o ^3;(t)/Aa;(0) ~ 
e'^*, where A is the Lyapunov exponent, Aa;(0) and Ax{t) are the uncertainties at times 
and t. When the Lyapunov exponent A < 0, ^(t) characterizes the rate of contraction 
towards periodic orbits. On the other hand, for A > 0, it characterizes the rate of divergence 
of chaotic orbits. At bifurcation and critical points (i.e., onset to chaos) the Lyapunov 
exponent A vanishes. Recently, it was shown that this feature is related to a power-law 
sensitivity to initial conditions on the form 

e(t) = [l + (l-g)A,t]V(i-^) , (1) 

with Ag defining a characteristic time scale after which the power-law behavior sets up. 

A quantitative way to measure the sensitivity to initial conditions is to follow, from 
a particular partition of the phase space, the temporal evolution of the number of cells 
W{t) occupied by an ensemble of identical copies of the system. For periodic and chaotic 
orbits W{t) = W{0)e'^^. In the particular case of equiprobability, the well-known Pesin 
equality reads K = XifX>0^ with K being the Kolmogorov- Sinai entropy |^ defined 
as the variation per unit time of the standard Boltzmann-Gibbs entropy S = — Z^Pjlnpj. 
This equality provides a link between the sensitivity to initial conditions and the dynamic 
evolution of the relevant entropy. 

At bifurcation and critical points and for an ensemble of initial conditions concentrated 
in a single cell, i.e. W{0) = 1, it has been shown that 

H^(t) = [l + (l-g)K,t]V(i-'?) , (2) 

with Kg being the generalized Kolmogorov-Sinai entropy defined as the rate of variation of 
the non-extensive Tsallis entropy Sg = {1 — J^pD/il ~ 1) 0- The Pesin equality can be 
generalized as Kg = Xg ii Xg > Tsallis entropies have been successfully applied to 
recent studies of a series of non-extensive systems and provided a theoretical background to 
the understanding of some of their unusual physical properties 

The expansion towards the critical attractor of an ensemble of initial conditions concen- 
trated around the inflexion point of the map can be characterized by a proper Sg evolving 
at a constant rate. Scaling arguments have shown that the appropriate entropic index q is 
related to the multifractal structure of the critical dynamical attractor by 

111 

(3) 



1 Q CXmin ^max 

where a^m and amax are the extremal singularity strengths of the multifractal spectrum of 
the critical attractor 0]. The above scaling relation has been shown to hold for the families 
of generalized logistic and circle maps P,p!oHT^. 
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However, the temporal evolution of critical dynamical systems can be strongly dependent 
on the particular initial ensemble. Although some scaling laws can be found for an ensemble 
of initial conditions concentrated around the map inflexion point, these are usually not 
universal with respect to a general ensemble. In this work, we are going to numerically 
investigate the critical temporal evolution of the volume of the phase space occupied by an 
ensemble of initial conditions spread over the entire phase space. This ensemble is expected 
to contract towards the critical attractor. Using a family of one-dimensional generalized 
logistic maps having df < 1, we will perform a detailed study of the parametric dependence 
of W{t) on the fractal dimension of the critical attractor. Due to the discrete scale invariance 
of the critical attractor, the convergence displays log-periodic oscillations [|13|. We are also 
going to explore the dependence of the amplitude of these oscillations with respect to the 
attractor 's fractal dimension. Further, the behavior of W{t) will be investigated for the 
one-dimensional critical circle map having df = 1. For this map, the temporal evolution is 
expected to display distinct trends since the critical attractor is dense. 



II. THE CONVERGENCE TO THE CRITICAL ATTRACTOR OF 
GENERALIZED LOGISTIC MAPS 

Logistic-like maps are the simplest one- dimensional nonlinear dynamical systems which 
allow a close investigation of a series of critical exponents related to the onset of chaotic 
orbits. This family reads 

xt+i = 1 - a\xt\' ; {z>l ; 0<a<2 ;t = 0,1,2,... ;xjG[-1,1]) . (4) 

Here z is the inflexion of the map in the neighborhood of the extremal point x = 0. These 
maps are well known to have topological properties not dependent of z. However, the 
metrical properties, such as Feigenbaum exponents fT^JT^ and the multifractal spectrum 
of the critical attractor do depend on z. In particular the fractal dimension of the critical 
attractor df{z) < 1 and therefore it does not fill a finite fraction of the phase space. 



For a set of initial conditions spread in the vicinity of the inflexion point, it was found that 
the volume in phase space occupied by the ensemble grows following a rich pattern with the 
upper bounds Wmax{t) governed by a power- law W^ax{t) oc t^/^^-''\ where q is the entropic 
index characterizing the relevant Tsallis entropy that grows at a constant rate. It has been 
shown that the dynamic exponent 1/(1 — g) is directly related to geometric scaling exponents 
related to the extremal sets of the dynamic attractor 0. 

Due to the presence of long-range spatial and temporal correlations at criticality, one ex- 
pects the critical exponent governing the temporal evolution to be sensitive to the particular 
initial ensemble. Indeed, the multifractal spectrum characterizing the critical dynamical at- 
tractor indicates that an infinite set of exponents are needed to fully characterize the scaling 
behavior. In particular, an ensemble consisting of a set of identical systems whose initial 
conditions is spread over the entire phase-space is a common one when studying non-linear 
as well as thermodynamical systems. 

Here, we will follow the dynamic evolution, in phase space, of an ensemble of initial 
conditions uniformly distributed over the phase-space and explore its relation with the gen- 
eralized fractal dimensions of the critical attractor. In practice, a partition of the phase 
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space on Ni,ox cells of equal size is performed and a set of Nc identical copies of the system 
is followed whose initial conditions are uniformly spread over the phase-space. The ratio 
r = Nc/Nj,ox is a control parameter giving the degree of sampling of the phase-space. 

Within the non-extensive Tsallis statistics, there is a proper entropy Sg evolving at a 
constant rate such that 



K,= \im [Sgit)-SM]/t (5) 

goes to a constant value as t —>■ oo. Notice that Kg < for the process of convergence 
towards the critical attractor. Assuming that all cells of the partition are occupied with 
equal probability, the entropy Sgit) can be written as 

q — I 1 — 

The last two equations imply that the number of occupied cells evolves in time as 

W{t) = [W{Oy~'^ + (1 - q)Kgtf'^-'''> (7) 

with the exponent /i = — — 1)>0 governing the asymptotic power-law decay. 

In figure 1, we show our results for W{t)/Nbox in the standard logistic map with inflexion 
z = 2 and from distinct partitions of the phase space with sampling ratio r = 0.1 . We 
observe that, after a short transient period when W{t) is nearly constant, a power-law 
contraction of the volume occupied by the ensemble sets up. W(t) saturates at a finite 
fraction corresponding to the phase space volume occupied by the critical attractor on a 
given finite partition. The saturation is postponed when a finer partition is used once the 
fraction occupied by the critical attractor vanishes in the limit Nhox oo. 

In figure 2, we show W{t)/Ni,ox for a given fine partition of the phase-space and distinct 
sampling ratios r. We notice that the crossover regime to the power- law scaling is quite 
short for large values of r so that a clear power-law scaling regime sets up even at early 
times. This feature is consistent with eq.(7) which states that the crossover time r scales 
as r ~ 1/1^(0)^-^ Further, the scaling regime exhibits log-periodic oscillations once the 
multifractal nature of the critical attractor is closely probed by such dense ensemble. A 
general form for W{t) reflecting the discrete scale invariance of the attractor can be written 
as 

wm^r^p[^) (8) 

where P is a function of period unity and A is the characteristic scaling factor between the 
periods of two consecutive oscillations. These log-periodic oscillations have been observed in 



a large number of systems exhibiting discrete scale invariance [|T3| . In general the amplitude 
of these oscillations ranges form 10~^ up to 10~^. Keeping only the first term in a Fourier 
series of P(lnt/lnA), one can write W{t) in the form 



W{t) = Cot" 



Ci / \nt 
2— cos 27r— + 
Co \ m A 



(9) 
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Log-periodic modulations correcting a pure power-law have been found in several sys- 
tems, as, for example, diffusion-limited-aggregation [|l^, crack growth |jl9||, earthquakes 
20| and financial markets It has also been observed in thermodynamic systems with 
fractal-like energy spectrum [^p3| . The factors controlling the log-periodic relative ampli- 
tude 2ci/co are not well known for most of the systems where it has been observed. In the 
present study, we can closely investigate the factors which may control these amplitudes by 
measuring it as a function of the map inflexion z for a fixed partition and sampling ratio (see 
figure 3). We found that these oscillations have amplitudes decaying exponentially with z 
as shown in figure 4. It is interesting to point out that the fractal dimension of the attractor 
is a monotonically decreasing function of z. Therefore, the above trend indicates a possible 
correlation between the amplitude of the log-periodic oscillations and the fractal dimension 
of the dynamical attractor. 

We also measured the critical exponent as a function of the map inflexion z. Our 
results are summarized in the Table. It is a decreasing function of z as can be seen in 
figure 5. The volume occupied by the ensemble depicts a fast contraction for 2; ~ 1 where 
the fractal dimension is small. On the other side, a very slow contraction is observed for large 
values of z, pointing towards a saturation or at most to a logarithmic decrease of Wif) in 
the limit of dense attractors. We would like to point out here that the exponent governing 
the expansion of the volume occupied by an ensemble of initial conditions concentrated 
around the inflexion point exhibits a reversed trend. Although scaling arguments have 
shown that this exponent can be written in terms of scaling exponents characterizing the 
extremal sets in the attractor, we could not devise a simple scaling relation between /i and 
the multifractal singularity spectrum. However, we observed that, when plotted against the 
fractal dimension of the attractor as shown in figure 6, the dynamic exponent [i is very 
well fitted by /i oc (1 — dfY ■, which indicates dj as the relevant geometric exponent coupled 
to the dynamics of the uniform ensemble. We would like to mention here that the same 
dynamic exponents were obtained for the generalized periodic maps which belong to the 



same universality class of logistic-like maps 110 



III. THE CONVERGENCE TO THE CRITICAL ATTRACTOR OF THE CIRCLE 

MAP 

The results from the previous section indicate that a slow convergence to the critical 
attractor shall be expected for dense critical attractors. However, it is not clear in what 
fashion this convergence will take place when the dynamical attractor fills the phase space 
with a multifractal probability density as occurs for the one- dimensional critical circle map 

e^^-^ = + _ ^ sin (27r^t) morf(l) (10) 

where < < 1 is a point on a circle. The circle map describes dynamical systems 
possessing a natural frequency lOx which are driven by an external force of frequency 102 
(f2 = UJ1/UJ2 is the bare winding number) and belongs to the same universality class of the 
forced Rayleigh-Benard convection [0. For Vt = 0.606661... the circle map has a cubic 



inflexion {z = 3) in the vicinity of the point ^ = 0. Starting from a given point on the circle. 
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it generates a quasi-periodic orbit which fills the phase-space and the dynamical attractor 
is a multifractal with fractal dimension df = 1 

In figure 7 we show our results for the temporal evolution of the phase-space volume 
occupied by an ensemble of initial conditions uniformly spread over the circle. W{t) exhibits 
a rich pattern which resembles the one observed for the sensitivity function associated to 
the expansion of the phase-space from initial conditions concentrated around the infiexion 
point. However, W{t) does not present any power-law regime. Instead, the lower bounds 
display a slow logarithmic decrease with time, saturating at a finite volume fraction. The 
saturation is a feature related to the finite partition used in the numerical calculation. 
This minimum decreases logarithmically with the number of cells in the phase-space as 
shown in figure 8. We also observed the same behavior for generalized circle maps with an 
arbitrary infiexion z |]12[. The critical attractor s within this family have aA\ df = 1 although 
they exhibit a ^-dependent multifractal singularity spectra. The ^-independent scenario for 
W{t) corroborates the conjecture that df is the relevant geometric exponent coupled to the 
dynamics of the uniform ensemble. 



IV. SUMMARY AND CONCLUSIONS 

In this work, we studied the temporal evolution in phase space of an ensemble of identical 
copies of one-dimensional nonlinear dissipative maps. We found that the phase-space volume 
occupied by an initially uniform ensemble displays a power-law decay with log-periodic 
oscillations whenever the dynamical attractor has a fractal dimension df < 1, i.e., when 
the fractal attractor does not densely fill the phase space. Generally, these oscillations 
also emerge in open high-dimensional systems operating at a self-organized critical state. 
The spatio-temporal long-range correlations present on the critical state refiect the scale 
invariance of the dynamical attractor. Therefore, the present work corroborates the concept 
that the fractal nature of the dynamical attractor and the presence of a characteristic scaling 



factor are key ingredients for the emergence of log-periodic oscillations 1 13]. The amplitude 
of the oscillations was found to depict a monotonic parametric dependence on df. For dense 
multifractal attractors, W{t) presents only a slow logarithmic contraction of its lower bounds 
followed by a rich pattern. 

The critical exponent characterizing the contraction of the uniform ensemble was found 
to have no direct relation to the one governing the expansion from a set of initial condi- 
tions concentrated around the inflexion point. In particular, no power-law was found for 
the contraction in the standard and generalized circle maps, in contrast to the ^-dependent 
power-law expansion. These results indicate that the relevant Tsallis entropy evolving at a 
constant rate (modulated by log-periodic oscillations) is characterized by an entropic index 
q that depends on the initial ensemble. It would be valuable to investigate the possible 
existence of classes of ensembles with a common dynamics in phase-space and, therefore, 
characterized by the same entropic index q. The non-universality of q with respect to the 
initial ensemble is related to the multifractal character of the dynamical attractor. However, 
as for the ensemble concentrated at the vicinity of the inflexion point, the exponent gov- 
erning the dynamics of the uniform ensemble is coupled to a geometric scaling exponent, in 
particular to the proper fractal dimension of the attractor. Extensive numerical work would 
be valuable to verify the validity of the proposed relation on higher- dimensional systems. In 
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any case, the present results come in favor of the concept that the degree of nonextcnsivity 
of the entropy measure evolving at a constant rate is related to the fractal nature of the 
dynamical attractor. 
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FIGURE AND TABLE CAPTIONS 



Figure 1 - The volume occupied by tlie ensemble W{t) (number of occupied boxes) as 
a function of discrete time in the standard logistic map {z = 2) and with sampling ratio 
r = 0.1. Prom top to bottom Nf,o^ = 2000, 8000, 32000, 128000. 

Figure 2 - The volume occupied by the ensemble W{t) as a function of discrete time in the 
standard logistic map {z = 2) and for a partition containing N^o^ = 128000 cells. Notice the 
emergence of log-periodic oscillations for large sampling ratios. 

Figure 3 - The periodic function W{t) / (cot""^) versus discrete time within the scaling regime 
and for r = 10. Data from map inflexions z = 1.1, 1.25, 1.5, 2.0 are shown. The amplitude of 
the oscillations decreases monotonically as z increases, but the characteristic scaling factor 
between the periods of two consecutive oscillations is roughly ^^-independent. 

Figure 4 - The amplitude of the log-periodic oscillations 2ci/co as a function of the map 
inflexion z for sampling ratio r = 10. The monotonic decrease of the oscillations indicates a 
close relation between these and the fractal dimension of the underlying dynamical attractor. 

Figure 5 - The dynamic exponent // governing the contraction of the occupied phase space 
volume [W{t) oc t~^^] as a function of the map inflexion z. 

Figure 6 - logj^o(A*) versus log^Q{df). The parametric dependence of the dynamic exponent 
II with the fractal dimension of the critical attractor is very well fitted to the form 
II oc {1 — dfY- It indicates that df is the relevant geometric exponent coupled to the 
dynamics of the uniform ensemble. 

Figure 7 - The volume occupied by the ensemble W{t) as a function of discrete time in the 
standard critical circle map. The lower bounds display a slow logarithmic decay with time 
saturating at a finite volume fraction due to the finite partition of the phase space. 

Figure 8 - The asymptotic lower bounds for the occupied volume in the phase space versus 
the number of cells Nbox- The logarithmic decay agrees with the prediction that ii{df 1) — > 
0. The same behavior was observed for the family of generalized circle maps and corroborates 
the conjecture that df is the relevant geometric exponent coupled to the dynamics of the 
uniform ensemble. 

Table - Numerical values, within the ^-generalized family of logistic maps, of: i) the 
dynamic exponent n governing the contraction towards the critical attractor of the uniform 
ensemble; ii) the entropic index q of the proper Tsallis entropy decreasing at a constant 
rate; iii) the fractal dimension df oi the critical attractor. These values also hold for the 
generalized periodic maps. The last line represents our results for the ^-generalized circle 
maps. 
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Table 
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1.10 
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1.62 ± 0.01 
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